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Abstract. Breather modes of the mKdV equation on the real line are known 
to be elastic under collisions with other breathers and solitons. This fact in- 
dicates very strong stability properties of breathers. In this note we describe 
a rigorous, mathematical proof of the stability of breathers under a class of 
small perturbations. Our proof involves the existence of a nonlinear equation 
satisfied by all breather profiles, and a new Lyapunov functional which con- 
trols the dynamics of small perturbations and instability modes. In order to 
construct such a functional, we work in a subspace of the energy one. However, 
our proof introduces new ideas in order to attack the corresponding stability 
problem in the energy space. Some remarks about the sine-Gordon case are 
also considered. 



1. Introduction 

Breathers have become a paradigm of beauty and complexity in nonlinear in- 
tegrable systems [T]. Although present in only some very particular models (sine- 
Gordon (SG), modified Korteweg-de Vries (mKdV), nonlinear Schrddinger (NLS), 
among others), their fascinating mixed behavior, combining oscillatory and soliton 
character, has attracted the attention of many researchers in the last thirty years 
[21 [3] . From the physical point of view, breather solutions seem to be relevant to 
localization-type phenomena in optics, condensed matter physics and biophysics 
[4]. They also play an important role in the modeling of freak and rogue waves 
events on surface gravity waves and also of internal waves in the stratified ocean, 
in Josephson junctions and even in nonlinear optics. See [5] for a representative set 
of these examples. 

Many authors have considered the study of dynamical properties of integrable 
systems. In addition to the existence of well-known soliton and multi-soliton solu- 
tions, and the elastic character of the corresponding interaction [I], one has to add 
that, from the Inverse Scattering Theory, the evolution of a rapidly decaying ini- 
tial data can be described by purely algebraic methods. Bounded energy solutions 
use to decompose into a superposition of radiation plus a very particular set of 
nonlinear elements (this is called the soliton resolution conjecture, see e.g. Schuur 
[2]). In such a decomposition, and for a very particular class of integrable models, 
not only solitons are allowed to appear, but also breathers and even more complex 
solutions. The emergence of solitons is not a surprise, since they are stable in the 
energy space [6], and even in less regular spaces, where only the mass is conserved 
0. 
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Stability of breathers 



However, in order to fully understand the aforementioned soliton resolution con- 
jecture, the stability and instability properties of breathers should be taken into 
account. Numerical computations suggest a positive answer to the question of the 
mKdV breather stability on the real line [8J. However, the simple problem of a rig- 
orous proof for their orbital stability has become a very challenging open problem. 

In this note, we announce that breathers solutions of the completely integrable 
mKdV equation [T| 

u t + (u xx + u 3 ) x = 0, (1.1) 

arc globally stable under suitable perturbations [9]. In more mathematical terms, 
we prove that breathers are orbitally stable for initial data which is a small pertur- 
bation of any breather solution, in a more regular subspace of the corresponding 
energy space. Our proof also works in the case of the sine-Gordon equation [1] 

utt — u xx + sin u = 0, (1.2) 

provided a suitable linearized problem has the adequate spectral properties. 



2. Breathers 



In what follows, every integral below is taken over R. 

First of all, and for the sake of completeness, we introduce the standard Sobolev 
spaces L 2 and H k : I? := {u : j u 2 < +00} , H 1 := {u £ L 2 : j u 2 x < +00} (the 
energy space), and H 2 := {u € H 1 : J u 2 x < +00} and so on. On the other hand, 
given a solution u(t) of (jl.l[) , the conservation laws at the H 1 -level of regularity 
are the mass 

M[u){t) :=\ J u* = M[u){Q), (2.1) 

and energy 

E[ U ]{t):=\Jul-\Ju i = E[um. (2.2) 

A satisfactory Cauchy theory is also present at that level of regularity, see e.g. 
Kenig-Ponce-Vega [TO]. Solutions u(t) constructed with initial data u(0) <G H 1 
are globally well-defined: one has u(t) G H 1 for all t, and mass and energy are 
conserved along the corresponding flow. 

Roughly speaking, breathers are time periodic, spatially localized solutions of 
(jl.lj) . Indeed, given a, j3 > 0, x\, x<z £ R, the breathers of mKdV (jl.lj) are given by 

„ , , T fP siri( a( x + St + xi)) \i . „. 

B a , p (t,x;x 1 ,x 2 ):=2V2d x arctan ^ V 1 lj , 2.3 

L \a cosh(p(a; + jt + Xi))J J 

with S := a 2 — 3/3 2 and 7 := 3a 2 — (3 2 (see [H]). Here a and (3 represent the two 
scaling parameters of the breather (note that a is also an oscillatory parameter) , 
and —7 describes the velocity of the solution. Note that B a ^ is periodic in time, 
but not in space, and this will be essential in our proof. If we take the limit /3/a « 1 
in (|2.3|) . this allows to simplify the expression for the breather to 

B ay p(t, x;0,0) ps 2V2/3cos(a(x + St)) sech((3(x +7*)) + o(-), 

and from a qualitative point of view, it shows explicitly its wave packet like nature, 
as an oscillation modulated by an exponentially decaying function. 
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We are ready to state the main result in [5] (here B a ^ refers to the breather 
solution 

Theorem 2.1 (instability of mKdV breathers). Let a, f3 > 0. There exist pa- 
rameters r]o,A , such that the following holds. Consider it G H 2 (R), and assume 
that there exists r\ G (0,770) such that 

||«o--Ba,/j(0,-;0,0)||H»(B) < V- (2-4) 

Then there exist Xi(t),X2{t) G IR such that the solution u(t) of the Cauchy problem 
for the mKdV equation U.l\) , with initial data uq, satisfies 

supllu(i) - B a ^(t, ■; x 1 (t),x 2 (t))\\ 2( < A r]. (2.5) 

tm y ' 

Finally, note that the constant Aq may depend on a and j3, but it is independent 
of n. There are explicit bounds on the variation of the parameters x\(t) and X2(t); 
see [5] for more details. 

3. Ideas of the proof 

In order to explain the main ideas involved in the proof of Theorem 12. 1[ let us 
start by recalling the proof for the case of mKdV solitons. First of all, solitons are 
regarded as minimizers of a constrained functional in the -ff^-topology. They are 
given by the expression 

u(t,x) = Q c (x-ct), Q c (s) ~ VcQ(Vcs), c>0, (3.1) 

with Q(s) := co ^ s - ) = 2v / 2<9 s [ai'ctan(e s )]. By replacing (|3.1[) in one has that 

Q c > satisfies the nonlinear ODE 

Q" -cQ c + Q 3 c = 0, Q c G H 1 . (3.2) 

A very useful method to guess either the stable or unstable character of a soliton 
is given by the following positivity criteria: since the mass of a soliton is given by 
M[Q C ] = then 

d c M[Q c ] = cT 1/2 > 0. 
Whenever the expression above is either zero or negative, the soliton is actually 
unstable [6j [12] . With this information in mind, one of the main ingredients of the 
arguments employed in some of the results mentioned above is the introduction of 
a suitable Lyapunov functional, invariant in time and such that the soliton is a 
corresponding extremal point. For the mKdV case, this functional is given by 

H[u](t) = E[u]{t)+cM[u](t), 

where c > is the scaling of the solitary wave, and E[u], M[u] are given in (|2.1[) - 
(12.21) . A simple computation shows that for any z(t) G H 1 small, 

H[Q C + z](t) = H[Q C }- J z(Q" c - cQ c + Q 3 C ) + Q[z](t) + 0(||«(t)||^i). (3.3) 

The zero order term above is independent of time, while the first order term is zero 
from p. 21) . It turns out that the second order term 

<z\m = \ J[zI + {c-zq 2 )z 2 ] 



4 Stability of breathers 

is positive definite modulo two direction^ related to the invariance of the equa- 
tion under shift and scaling transformations. Time-dependent parameters are then 
introduced in order to remove these unstable modes. Once these directions are 
controlled, the stability property follows from (|3.3|) . 

In [9] , we follow similar lines to prove the instability theorem of mKdV breathers 
12.11 First of all, recall the mass and energy of breathers: 

M[B a , p ] = 4/3, E[B a , p ] = (3.4) 

From the first identity, standard stability tests can be deduced: 

d a M[B a ,p] = 0, dpM[B ai p] = 4 > 0, 

which suggest a kind of critical character of the scaling parameter a, and the 
subcritical behavior of /3, as expected. However, we will show in the next lines 
that, with enough regularity on hand, a behaves as a stable direction. 

The first step into that direction is the following. We have found an explicit, 
nonlinear elliptic equation satisfied by any breather B = B at p, for any fixed t £ K, 
which is the following: 

B (4x) - 2(/? 2 - a 2 )(B xx + B 3 ) + (a 2 + p 2 ) 2 B + 5BB 2 X + 5B 2 B XX + ^B 5 = 0. (3.5) 

The proof of this result is involved and requires several new identities. In particular, 
a key step for the proof is to establish the simpler second order identity 

B xt + 2B f BB t = (a 2 +[3 2 ) 2 B + 2((3 2 -a 2 ) f B t . 

J — oo J — oo 

It appears that (|3 . 5[) cannot be directly obtained from the original ideas of Lax [13] , 
since breathers do not decouple into solitons as time evolves. 

With the identity (|3 . 5|) on hand, one can compute explicit expressions for the 
linear operator associated to it. Indeed, if B\ :— d Xl B a ^ and B-i := d X2 B Ut p, 

C[B t ] = C[B 2 ] = 0, 

where C is the unbounded, self-adjoint operator 

C[z](x; t) := z ( 4*)(z) - 2(/3 2 - a 2 )z xx {x) + {a 2 + (3 2 ) 2 z{x) + 5B 2 z xx {x) 

+10BB x z x (x) + [5B 2 X + 10BB XX + yB 4 - 6(/3 2 - a 2 )B 2 ]z(x), 

(3.6) 

defined in the dense subspace H 4 of L 2 . Additionally, if A a B := d a B and ApB := 
dpB, one has 

C[A a B] = -4a[B xx + B 3 + {a 2 + f3 2 )B], (3.7) 
C[A a B] = 4f3[B xx + B 3 - (a 2 + (3 2 )B}. (3.8) 

Finally, if we define B := "stff^+^f ' we § et 

C[B ] = -B. (3.9) 
These expressions will be useful in last part of the proof. 



"^A quadratic form Q has a negative direction z ^ if Q[z] < 0. Note that Xz is also a negative 
direction, for all real-valued A 7^ 0. 
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A detailed study of the linear operator C (j9|) reveals that, for any t £ K its 
kernel is isolated and it is spanned by B\ and B2, it has continuous spectrum given 
by the interval [(a 2 + /3 2 ) 2 ,+oo) in the case j3 > a, and [4a 2 /Si 2 , +00) in the case 
P < a. Even more surprising, it has just one negative eigenvalue. The first hint 
of this result can be obtained by computing the quadratic form associated to C for 
the elements which are the first candidates to be negative directions, namely the- 
scaling parameters A a B and ApB. Indeed, from (|3.7j) and (|3.4j) one has 

J C[A a B]A a B = Aad a E[B] = 32a 2 /3 > 0, 

J £[Af}B]A B = -4:/3df 3 E[B} - 16/3(a 2 + ft 2 ) = -32a 2 [3 < 0. 

A rigorous proof of the existence of a unique negative eigenvalue makes use of 
the spectral theory developed by L. Greenberg [14]. As a conclusion, it is expected 
that, as in the standard soliton case, the scaling parameter j3 describes an instability 
mode, which need to be removed, or well-estimated. 

Finally, the more important ingredient appearing in the proof of 12.11 in [9] is 
the introduction of a Lyapunov functional, well-defined in the H 2 topology, and for 
which breathers are surprisingly not only extremal points, but also local minimizers, 
up to symmetries. This functional is a suitable combination of the mass and energy 
(|2.ip - ([2.2[) . and a third conserved quantity, defined at the 7J 2 -level by 

F[u](t) :=lju 2 xx -y u 2 u 2 x + \Ju & . (3.10) 

(The Cauchy theory for H 2 initial data is well understood and F remains constant 
along the flow.) We define this new Lyapunov functional as follows. Let B = B a .p 
be any breather with scalings a and /3, and M[u] and E[u] introduced in (|2.1[) . 
(|2.2|) . We denote, for any time t, 

H[u] := F[u] + 2(f3 2 - a 2 )E[u] + (a 2 + /3 2 ) 2 M[u}. (3.11) 

This functional is reminiscent of the one appearing in the foundational paper by 
Lax [13] , concerning the 2-soliton solution of the Korteweg-de Vries (KdV) equation, 
and generalized to the KdV AT-soliton states by Maddocks-Sachs \§\ • This idea has 
been successfully applied to several 2-soliton problems, for which the dynamics 
decouples into well-separated solitons as time evolves, see e.g. [15]. However, there 
was no evidence that this technique could be generalized to the case of even more 
complex solutions, such as breathers. Compared with those results, our proofs are 
involved, since there is no mass splitting as t — > +00. 

Coming back to the proof, it is clear that H[u] is a conserved quantity, well- 
defined for £f 2 -solutions of (|1.1|) . Moreover, one has the following decomposition: 
let z £ H 2 be any function with sufficiently small norm, and B = B a _p be any 
breather solution. Then, for all t £ M, one has 

H[B + z]-H[B] = ^Q[z]+Jf[z], (3.12) 

with Q being the quadratic form defined by L, namely Q[z] := J z£[z], and Af[z] 
satisfies |7V"[z]| < ||^|||/2- Note that in this result, the linear term in z vanishes since 
it is proportional to the left hand side of (|3.5|) . 
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Stability of breathers 



From the previous analysis, we see that the functional p. lip will allow us to 
control the dynamics of the perturbative term z(t), provided we manage the three 
instability directions that appear as a consequence of the symmetries satisfied by 
(|1.1[) : the two dimensional kernel Bi,B2, and the eigenfunction corresponding to 
the negative eigenvalue of C. We vary in time the parameters x\ and xi in (|2.3I) to 
satisfy, for all time, the orthogonality conditions 

J Biz = J B 2 z = 0, (3.13) 

with z(t) := u(t) — B a ^{t;xi(t),X2{t)). This is an absolutely necessary condition 
in order to obtain an orbital stability property. However, we do not modulate 
the scaling instability ft. Instead, we control the dynamics by first replacing the 
corresponding negative mode by a more tractable direction, the breather itself (this 
technique was first introduced by Weinstein in |16j): one has 

Wmfm + \\4t)\\ 3 H z>Qim>Mt)\\h- 

provided (|3.13[) is satisfied. Roughly speaking, this property is consequence of ()3.9j) 
and the fact that J BqB ^ 0, that is, Bq and B are not orthogonal, and Bq is 
also a negative direction. Finally, we use the mass conservation law to give a new 
estimate on the last term in (|3.14[) . which is at the first order independent of time, 
and only depending on z(0). This last element of the proof allows to control the 
term z(t) for all times, proving the stability property. 

4. Discussion 

Our results emphasize, and demonstrate, at the rigorous level, some deep con- 
nections between breathers and the 2-soliton solution of mKdV. This connection 
comes from the structure of the Lyapunov functional p. lip that we have used as 
an essential tool in the proof. The main theorem is stated for an H 2 perturbation, 
although we believe that it can be improved to reach the H 1 level, but with a harder 
proof. We also believe that our proof gives some insights to approach the study of 
the instability of periodic arrays of breathers, as numerically stated in |17 [ 118 ) . In 
order to prove this fact, it is necessary to find a suitable Lyapunov functional for 
which lattice breathers are critical points. 

Concerning the SG case, our proof also applies with no significant modifications 
(treating SG as a matrix operator problem for [u, ut) £ H 2 x H ), provided a suit- 
able linearized operator (of matrix type and of fourth order in terms of derivatives), 
has only one negative eigenvalue, and the associated kernel is nondegenerate. For 
the moment, a proof of that result has escaped to us. Therefore, if these conditions 
are satisfied, SG breathers are H 2 x H 1 stable. 

However, it is highly expected that breathers could not survive under nontrivial 
perturbations of the equation, as is showed in [3| , in the SG case. 

5. Conclusion 

In this note we have presented some details of the proof of the nonlinear stability 
of breathers for the mKdV equation. Our conclusion is that mKdV breathers are 
nonlinear stable at the H 2 level of regularity. No systematic work on this kind was 
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reported before, probably as a direct consequence of the difficulties that arise by the 
complexity of the breather solution and the computations involved. As we explained 
above, the introduction of the Lyapunov functional f|3 . 1 f [) for which breathers are 
local minimizcrs, up to symmetries, is the essential point. This functional allows to 
control the perturbative terms and the instability directions that appear during of 
the dynamics. Finally, we indicate that our arguments are general and, in principle, 
can be applied to several equations with breather solutions. We expect to consider 
some of these problems in a forthcoming publication. 
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